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Abstract: Contour integrals of rational functions over A4o,n, the moduli space of 
n-punctured spheres, have recently appeared at the core of the tree-level S-matrix of 
massless particles in arbitrary dimensions. The contour is determined by the critical 
points of a certain Morse function on M. Q n . The integrand is a general rational function 
of the puncture locations with poles of arbitrary order as two punctures coincide. In 
this note we provide an algorithm for the analytic computation of any such integral. 
The algorithm uses three ingredients: an operation we call general KLT, Petersen’s 
theorem applied to the existence of a 2-factor in any 4-regular graph and Hamiltonian 
decompositions of certain 4-regular graphs. The procedure is iterative and reduces 
the computation of a general integral to that of simple building blocks. These are 
integrals which compute double-color-ordered partial amplitudes in a bi-adjoint cubic 
scalar theory. 
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1 Introduction 


The complete tree-level S-matrix of a large variety of field theories of massless particles 
are now known (or conjectured) to have a description in terms of contour integrals 
over M 0 ,n, the moduli space of n-punctured Riemann sphere [1-13]. Some of these 
theories are Yang-Mills, Einstein gravity, Dirac-Born-Infcld, and the U(N) non-linear 
sigma model [12, 13]. The new formulas for the scattering of n particles are given as a 
sum over multidimensional residues [14] on J\Ao, n . 

The position of n punctures on a sphere can be given using inhomogenous coordi¬ 
nates as {<7i, cr 2 ,..., a n }. Three of them can be fixed using PSL( 2, C) transformations, 
say fly, 02 , er 3 . Therefore the space is n — 3 dimensional and we are working locally on 
a patch isomorphic to C n_3 . The next step in the construction is a rational map from 
C n ~ 3 —» C n ~ 3 which is a function of the entries of a symmetric n x n matrix, s a b, with 
vanishing diagonal, i.e., s aa = 0, and all rows adding up to zero. These are the coordi¬ 
nates of the space of kinematic invariant for the scattering of n massless particles. The 
explicit form of the map is (oq, cr 5 , ..., a n } — » {E 4 , E 5 ,..., E n } with 


E a (a) 


£ 

b=l,by^a 


Sab 

&a &b 


for a G { 1 , 2 ,..., n}. 


( 1 . 1 ) 


Using this map, scattering amplitudes, denoted as M n , are defined as the sum over 
the residues of 

,2 H(a,k,e) 


]Qd(T a |123|- 


a=4 


E 4 (a)E 5 (a) ■ ■ ■ E n (a) 


( 1 . 2 ) 


over all the zeroes of the map {E 4 , E 5 , ..., E n }. Here 1123| = (cxi — a 2 ) (o ~2 ~ C 3 ) (03 — 0"i) 
and H(a, k, e) is a rational function that depends on the theory under consideration and 
contains all information regarding wave functions of the particles such as polarization 
vectors and momenta k The equations defining the zeroes, E 4 = E 5 = ■ ■ ■ E n = 0, 
are known as the scattering equations [15-23] More explicitly, 


M n = 


|123| 2 H(a*,k,e) 


, , f d(E 4 ,...E n ) 
a*ez(E) det ^ dia4 ^ an) 


(1.3) 


where Z(E) is the set of all zeroes of the map. This representation of scattering 
amplitudes is known as the Cachazo-He-Yuan (CHY) approach [1, 2, 11, 12] 
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The zeroes are generically isolated and are the values of a's for which the Morse 
function on A4 0 , n 

= ^ s ab In Wa ~ cr b | 2 (1.4) 

a<b 


has local extremes 1 [13, 18]. 

In this paper we are not concerned with particular theories. Instead, our aim is to 
provide an algorithm for the analytic computation of any integral of the form 


n da ° 

a =4 


11231 


E 4 (a)E 5 (a) ■ ■ ■ E n (a) 


F(a), 


(1.5) 


where T is the same contour as above, i.e., a sum over all residues at Z(E). Here F(<x) 
is any rational function of only the puncture coordinates cr’s which transforms as 

n . 

F{&) -A TT(c(T a + d) 4 F(<r), under a a -> , (1.6) 

t Li ca a + d 


with ad — be = 1, i.e., under an PSL( 2, C) transformation. 

The transformation of F(cr) ensures that the integral (1.5) is independent of both 
the choice of which puncture coordinates to fix and their values. The transformation 
also implies that F(a) is only a function of differences a a — <7 b which we denote as a ab - 
Clearly a ab = — a ba . The only other condition we impose on F(cr ab ) is that all its poles 
are of the form for some integer m > 0. 

The simplest kind of integrals are defined in terms of the so-called Parke-Taylor 
factors [25] defined for a particular ordering of n labels (a(l)a(2) • • • a(n)) with a G S n 
as 


(1.7) 


(a(l)a(2) • • • Q)(/l)) 0"a(l)o(2) &a(2)a(3) ' ' ' &a(n— l)a(n) < ^"a(n)ci( 1) 

Clearly, any Parke-Taylor factor has half the PSL( 2, C) weight needed to construct a 
valid F(a ab ) [1, 12], One can define integrals labeled by a pair a permutations a, (3 E S n 
using 


F a ’P(a ab ) = 


x 


a(l)a(2) ■ ■ • a(n)) (/3(1)P(2) ■ ■ ■ P(n)) ’ 


( 1 . 8 ) 


1 A Morse function is a real function with non-degenerate critical points [24]. 
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or more explicitly [2] 


t \p\ _ f j 1 1 

m(aW “ Jr ^ (a(l)«(2) •••«(«)) (/3(l)/3(2) • • ■ fS(n))' 

where we have introduced a shorthand notation for the measure 

11231 2 


(1.9) 


d/ y — d(J Q 


a=4 


E 4 (a)E 5 (a) ■ --Enja)' 


( 1 . 10 ) 


Integrals of the form m(a|/3) have been studied in the literature and are known 
to evaluate to a sum over connected tree Feynman graphs with only cubic (trivalent) 
interactions which are compatible with the two planar orderings defined by a and (3 
[2], We review this result in detail in section 2 and explain how to explicitly evaluate 
them as a rational function of the variables s a b- Here it suffices to say that these known 
integrals form the basic building blocks of our construction and the main result of this 
work is an algorithm for writing 


dn n F(cr ab ) = R(m(a\(3)), 


( 1 . 11 ) 


where R is a rational function of its variables with only numerical coefficients. 

The reason general integrals are of interest can be seen, for example, in the evalu¬ 
ation of an n graviton amplitude which contains a term of the form [1, 2, 11] 



(6i-e 2 ) 2 

a l2 


n / n 

n e 

a =3 \b=2,b^=a 


€a ■ h, 


&lb 

@ab@la 


2 


( 1 . 12 ) 


In this formula e c , k c are fixed data and after fully expanding (1.12) they can be factored 
out leaving arbitrarily complicated integrals of the form (1.11) to be evaluated. Also 
motivated by the same physical problem, Kalousios developed a technique, different 
from the one presented here, for the computation of general five-point integrals in [26]. 

The algorithm we develop is based on three key constructions. The first is a 
generalization of the Kawai-Lewellel-Tye (KLT) relation [27-29]. The KLT relation 
was originally discovered as a relation among closed and open string theory amplitudes 
but since then it has inspired similar relations in held theory and more recently it 
found a natural set up which allows vast generalizations in the CHY representation of 
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amplitudes. We present the general KLT construction in section 3. 

The second result is a classic one from graph theory [30]. Consider an integrand 
F(cr) such that it does not have any zeroes. This means that it is only the product of 
2 n factors a a b in the denominator with a trivial numerator that can be set to unity. 
Representing each puncture by a vertex and each a a b by an undirected edge connecting 
vertices a and b one finds that each F(a) leads to a unique 4-regular graph Gf (not 
necessarily simple). A classic result of Petersen guarantees that any 4-regular graph 
with n vertices is 2-factorable. This means that Gp it is always the union of two 
2-regular graphs with n vertices. Petersen’s result is reviewed in section 4. 

The third and final ingredient is an observation regarding the existence of a Hamil¬ 
tonian decomposition of graphs [30, 31]. In order to state the observation let us choose 
any 2-regular multigraph 2 G with n-vertices and no loops. We say that a connected 
2-regular graph with n-vertices, H conn , is compatible with G if the 4-regular graph 
obtained from the union of G and H conn contains two edge-disjoint Hamilton cycles. 
The observation is that out of the (n — 1)! possible connected graphs the number of 
compatible graphs with G is always larger than (n — 3)!. This is explained in section 5. 

In section 6 all ingredients are combined to produce the final algorithm for com¬ 
puting the rational function R in (1.11). The algorithm is general but in particular 
cases it can be modified to make it much more efficient. 

Section 7 is devote to examples that not only illustrate the use of the algorithm 
but also give the explicit Hamiltonian decompositions needed for the computation of 
the most general six-point integral. 

In section 8 we end with discussions including future directions and physical appli¬ 
cations in the form of novel relations among amplitudes. The appendix has a detailed 
explanation of how to implement Petersen’s theorem. The implementation is not far 
from being the actual proof so it a good way to gain intuition on why the theorem 
holds. 

2 In this work we use the terminology graph and multigraph interchangeably. In fact, the restriction 
to simple graphs is never necessary. 
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2 Definition of Building Blocks 


The aim of this work is to provide an algorithm for the reduction of contour integrals 
on the moduli space of an n-punctured sphere of the form 


J dfi n F(a ) 


( 2 . 1 ) 


in terms of a basis of known integrals. Ensuring that the integrand is PSL( 2, C) 
invariant implies that F(a) has the form 

F[a) = (12 ■ • • n )(T(l)^(2) • - • T(rz)) (2 ' 2) 

where (12 • • • n) is the canonical Parke-Taylor and ( 7 (l)y( 2 ) • • • y(n)) is a Parke-Taylor 
factor with a 7 G S n ordering (see (1.7) for the Parke-Taylor factor definition). / is a 
rational function of r^i which are general cross ratios, i.e., 


Tijkl — 


& ij ® kl 
Oil®jk 


(2.3) 


Of course, the choice of Parte-Taylor factor is completely arbitrary and can be conve¬ 
niently made depending on the case. The measure d/i n was defined in (1.10) and is 
reviewed below. 

In this section we discuss the basic building blocks which are special contour inte¬ 
grals with f{rijki ) = 1 and whose values are explicitly known [2, 4, 9]. The building 
blocks are labeled by a pair of permutations a, /3 E S n /Z n . The reason one has to mod 
out by cyclic permutations Z n is obvious from the definition 


m(a\j5) 



1 1 

(a(l)a(2) ■ ■ ■ a(n)) (/3(l)/3(2) ■ ■ • /3 (b)) ' 


(2.4) 


Recall that (ck( 1 )ck( 2 ) • • • ck(ti)) — < ^"a(i)a(2) ( ^"a(2)a(3)''' &a(n—i)a(n)&a(n)a(i) and the mea 
sure is 




11231 


a =4 


E 4 E 5 ■ ■ ■ E n 


(2.5) 


An explicit evaluation of the integral m(cx|/3) would involve solving the equations [4, 8 , 
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9, 26, 32, 33] 


E a (cr) — Y ——— = 0 for a 6 {4, 5,..., n}. (2.6) 

, °a — &b 

b=l, 

These equations have (n — 3)! solutions as proven in [1, 32] and the data s a b can be 
taken to be the components of a symmetric n x n matrix of complex entries such that 
Sll — $22 — ' ' ' s n n — 0 an d 


y S ab = 0 for a G {1, 2,..., n}. (2.7) 

b=l,b^a 

Once the solutions are found one computes the Jacobian matrix 




Sab 


E n 

c=l,c^a 


a 7 ^ b, 
a = b. 


( 2 . 8 ) 


Defining $(23 as th e ( n — 3) x (n — 3) minor of <f> obtained by deleting rows and columns 
1, 2, 3, one has that 


m(a\/3) 


(n—3)! 

E 


1=1 


1 123 1 2 1 1 

det (a(l)a( 2 ) • • • a(n)) (fi(l)fi( 2 ) ■ ■ -fi(n)) 


(2.9) 


where aa' 1 denotes the value of cr a on the I th solution. 

Solving the equations E a = 0 is a nontrivial task when n > 5 as for generic values 
of s a b an d after finding a Groebner basis one is faced with an irreducible polynomial of 
degree (n — 3)!. 

Luckily, it is easy to make a simple proposal for what m(a\{3) evaluates to and then 
prove that it is the right answer. This was done in a series of papers [2, 11, 32], Here 
we simply quote the result and use these integrals as building blocks for generic ones. 

2.1 Evaluating m(a\/3) 

Consider any connected tree graph T with n vertices of degree 3 one and n —2 vertices of 
degree three. Associating a label { 1 , 2 ,..., n} to the vertices of degree 1 one can assign 

3 The degree of a vertex is defined as the number of edges incident to the vertex. 
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a rational function of s a b to T as follows. Every internal edge e, i.e. not connected to 
a degree one vertex, divides the graph T into to graphs T L and T R if the edge e was 
removed. Let the subset of vertices of degree one from {1,2,..., n} which he on Tl be 
Sl and those on T R be S R . Then it is easy to show that 

Sab= Sab (2.10) 

a,fesS_L a,b£Su 

as a consequence of (2.7). Therefore this is a quantity that can be associated with the 
edge e and we denote it as 

P ,? = £ (2.H) 

a,6eS_L 

The reason for the notation is that in physical applications this is the norm of a Lorentz 
vector. 

The rational function associated with the graph T is then 

( 2 . 12 ) 

eG4"‘ 6 

where E™ 1 is the set of all internal edges of T. In physics terminology, T is a Feynman 
diagram in a massless cubic scalar theory and w(T) is the value of the graph obtained 
by using Feynman rules. 

A given diagram T can be drawn on a plane in a variety of ways. Each way of 
doing so defines a cyclic ordering of the labels {1,2 ,, n}. We say that T is consistent 
with an ordering a € S n /7* n if a is one of the possible orderings obtained when T is 
drawn on a plane. Let us denote the set of all graphs T consistent with the ordering a 
by T(a). 

Now we can state the main result of this subsection. The integral 

m(a\j3) = (~l) q Y W ( T )> ( 2 - 13 ) 

Ter(aOnr(/?) 

where q was defined in [2] and will not be relevant for our purposes 4 . 


4 For more details see the equation (3.4) in [2]. 



2.2 Examples 

Let us consider some simple examples in order to illustrate the use of the general 
formula (2.13). The first one is the four-point integral with the canonical ordering 

m(1234 | 1234) = — + —, (2.14) 

Sl2 S14 


which is a trivial computation using (2.12) and (2.13). One can also find different 
orderings a and /3 such that its result is just one Feynman diagram, for example 

m(1234| 1243) = —, (2.15) 

Sl2 


which is the diagram 



(Fig.2.1) Feynman diagram contributing for m(1234| 1243). 


More interesting examples are the five-point computations. For instance, with a 
and f3 in the canonical ordering one obtains 

m( 12345 112345) = —(2.16) 

S12S45 S12S34 S23S15 S23 s 45 S15S34 

In hve points we can also have two different orderings with intersection on only one 
Feynman diagram, for example 

m(12345 112534) =-, (2.17) 

S12S34 


with diagram 




3 


4 


(Fig.2.2) Feynman diagram contributing for m(12345 | 12534). 
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In addition, the m(a\/3) five-point integrals have more properties than the m(a\/3) four 
point matrix. For example, one can think in a Parke-Taylor with an a ordering as 
a vector and the matrix element m(a\/3) as the inner product among two of them. 
So, a natural and interesting question arises, given a Parke-Taylor with a particular 
ordering, what is its orthogonal space? for instance, at four points, there does not exist 
two orthogonal orderings. But, in five points, every Parke-Taylor with an a ordering has 
a 1-dimensional orthogonal space. For example, let us consider the canonical ordering 
(12345), its orthogonal space is generated by the Park-Taylor (14253), i.e 

m(12345 114253) = 0. (2.18) 


3 Generalized KLT 

In this section we introduce the first result needed for the computation of general 
integrals. In the 80’s Kawai, Lewellcn and Tye (KLT) found a relation connecting 
scattering amplitudes of closed strings to the sum of products of open strings amplitudes 
[27-29]. While closed string amplitudes are computed on the sphere and hence are 
permutation invariant, open strings are defined as sums over partial amplitudes. Each 
partial amplitude is computed on a disk where the external states are inserted on the 
boundary and therefore possess an ordering. Thanks to the Bern-Carrasco-Johansson 
relations (BCJ) [34], the modern version of the KLT formula can be written as (for 
more details see [35-39]) 

M do„d = V (3.1) 

a,/3eS„_ 3 

where permutations of {2, 3,..., n — 2} are denoted a and j3. In the formula above, 
stnng( q, |y 3 ) j s called the KLT momentum kernel [40] and it is a somewhat complicated 
function of the variables s a b whose explicit form is not relevant at this point. 

The string theory formula (3.1) has a field theoretic analog obtained by taking the 
infinite tension limit and it relates amplitudes of gravitons to that of gluons. More 
explicitly [28, 29, 34, 38, 40], 

Mr“= £ M5 1 "“(l,Q,n, n -l)5(<S|/5)M?™(l,ftn-l,n). (3.2) 

a , 0 £ Sn -3 
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Here <S(d|/3) is the infinite tension limit of 5 strmg (d|/3). 

In [12], it was realized that the field theory KLT relation (3.2) and many general¬ 
izations naturally follow from the CHY representation of amplitudes. 

Let us summarize the construction with special emphasis on the structures needed 
in section 6. 

Consider any contour integral as a starting point (sp), it plays the role of M n in 
(3.1), 

h P = j dfi n :1(a), (3.3) 

with an integrand that can be separated into two parts 


I{<t)=I l {<t)I r {<t), (3.4) 

where each “half-integrand” has half the PSL( 2, C) weight of the full integrand. Ex¬ 
amples of such integrands and half-integrands were studied in the previous section with 
Parke-Taylor factors being the half-integrands. 

The evaluation of (3.3) is given by 

(n_3)! I123I 2 

Let us denote the combination det ^ill/l 123| 2 as det 7< L. 

Introducing an (n— 3)!x(n—3)! diagonal matrix in solution space Du = det 7 ^(aa^)Su 
and (n — 3)!-dimensional vectors II and Tr one finds a matrix form of (3.5) 

L sp = Xl D 1 Xr- (3.6) 


T (P 


(3.5) 


The next step is to find an alternative representation of the matrix D in terms of a 
(n — 3)! x (n — 3)! matrix but this time in the ordering space. A natural candidate is 
to consider a submatrix of the (n — 1)! x (n — 1)! matrix whose entries are given by 
m(a\/3) with a,/3 E S n /Z n . 

More explicitly, the definition in (2.4) gives 


m(a |/3) 


(n—3)! 

E 


7=1 


11 1 
det 7 $ (a(l)a(2) ■ ■ • a(n)) (/3(l)/3(2) ■ • • /3(n)) 


(T a =(T 


{!) ' 


(3.7) 


11 







This time it is convenient to introduce a rectangular (n — 1)! x (n — 3)! matrix Q with 
entries 

(3.8) 


^ a (a(l)a(2) • • • a(n)) det' $ 


x 


M) 


so that 


(n—3)! 


™(a |/3) = Qi D uQ 


(3.9) 


i,j =l 


Next we divide the discussion into the derivation of the standard KLT result (3.2) 
and then its most general form. 

3.1 Standard KLT 

Let us first discuss how to recover the standard KLT formula in its modern version 
[27-29, 34, 38] before proceeding to the more general discussion. As discussed above 
one has to select an (n — 3)! x (n — 3)! submatrix of the matrix m(a |/3). The choice 
that leads to the modern version of the KLT formula is 


m(l, a, n — 1, n|l, (3, n, n — 1) = (d|/3), 


(3.10) 


with a, (3 G S n -3 permutations of the remaining n — 3 labels. Once the choice has been 
made two square matrices can be defined 


Ul = Ql A ,n,n- X, V^Q[ A n-l, 


l — nf 


(3.11) 


Finally, it is possible find a representation for D 


m KLT = U l 'DV D = (V) l m KLI ((; 1 ) 




(3.12) 


Using this formula in (3.6) 


J sp = 1 L D~ l X R = (UX L y ( m ^ 1 )' 1 VI 2 


R- 


(3.13) 
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It is easy to recognize that 


(Ul L )(l,a,n, n-l) = [ dp n Xl ^ , (3.14) 

J (1, a,n,n — 1) 

(yi R )(lJ,n-l,n)= [ rf/v 

J (1, p, n — 1, n) 


When applied to gravity and Yang-Mills amplitudes (3.13) becomes the standard KLT 
formula (3.2). Even though it is not used in this work, let us make this more explicit 
for completeness. In [1] gravity and Yang-Mills amplitudes are computed as follows 


Mf avitons = J ^ 

M| luon » = [ dp 


(Pf / 'h(e, k, cr)) 2 , 

Pf /v k(e, k, a) 
(a(l)a(2) • • • a(n)) ’ 


(3.15) 


where 'I'(e, k , cr) is some 2 n x 2 n matrix whose precise form can be found in [1], Now 
it is clear how this leads to the KLT formula directly. Moreover, it shows that the 
momentum kernel S(a\f3) = (m KLT ) -1 (d|/3). 


3.2 General KLT 

Let us now discuss the more general construction. It is clear that the same steps can be 
followed as in the standard KLT construction if one chooses a general (n — 3)! x (n — 3)! 
sub-matrix of the matrix m(a\/3). However, not all sub-matrices are allowed as the 
construction requires the computation of its inverse. It turns out that m(a\/3) has 
vanishing determinant and so do some of its (n — 3)! x (n — 3)! sub-matrices. 

Let C and 1Z be both subsets of permutations S n /Z n with (n — 3)! elements. We 
say that C and 1Z are independent if the matrix with entries 


m c \ n = {m(a\/3) : a E C, /3 E 1Z} 


(3.16) 


has non-vanishing determinant. 

Provided £ and 1Z are independent one can obtained a formula for / sp of the form 

I. p= E (™ £ ’ K )i (VI K )ft (3.17) 
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where we have defined 


(3.18) 


U I a = Q I a > a e C, 
V^Q 1 ^ 6e U, 


and 


m c ) a 

<yi K ) 3 



IlM 

(a(l)a( 2 ) • • • a(n)) 

IrW) 

(pmm-m)’ 


cx G 

pen. 


(3.19) 


This is the most general form of the KLT relation that is needed in the algorithm 
presented in section 6 . 


3.3 Examples 

In this subsection we give a simple example to show how the standard KLT construction 
can be used in the computation of residue integrals. After that, we formulate an 
example where the standard KLT is not enough. 

Let us start with an example where the standard KLT construction suffices. Con¬ 
sider the five point integral 

f 1 i 

/sp = ,/ (12345) (12) (345)' (3 ' 20) 

Recall the general definition given in the introduction which applied to this case implies 
(12) = <Ti 2 cr 2 i and (345) = < 734045 ( 753 . This integral is not of the form studied in section 
2. The idea is then to fold a way of writing it in terms of the building blocks of section 
2 . 

Let the KLT basis be given by permutations (l,b,32) and (l,d, 23) with a per¬ 
mutations of {4,5}. 

Using the KLT formula one has 


/ sp = ^ m(12345|l,d,32)(m KLT )-^ 

Q!j/3Eperm(4,5) 


dn 5 


1 

(12)(345)(1, /3,23)' 


(3.21) 


At first sight it seems that the problem has been made worse as one has to now deal 
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with two new integrals 


J dfi5 ( 12 ) (345) (14523) ’ 

However, it is simple to check that 
(12) (345) (14523) = (12354) (12543), 


/ rf/i5 (12) (345) (15423)' 

(12) (345) (15423) = (13542)(12345). 


(3.22) 


(3.23) 


This means that we have succeeded in expressing I sp in terms of the building blocks. 
Using the explicit form of the building blocks one finds that 


dfi ,5 


(12345) (12)(345) 


E (&£)< 


( mKlT bi 


(3.24) 


a,/ 3 Sperm( 4 , 5 ) 

1 1 
- + - 


+ 


S 15 


+ 


Sl4 


+ 


Sl 5 


S12S34 S12S45 S? 2 S34 sf 2 S45 S12 s 45 


where the m KLT (l, a, 32|l,/3, 23) matrix is given by 


m KLT (l 

and the vectors 
(UX L )(1, a, 32) = 
(Vl R )(lJ, 23 ) = 



f 1 1 

1 \ 

11,/3,23) = | 

[ S23 S14 523 545 

1 1 

523 545 I 

1 1 / 


\ 523 545 

S 23 $15 $23 S45 / 


(3.25) 


11111 
-+-+-+-+- 


S12S45 S23S45’ S12S45 S12S34 S23S15 S23S45 S15S34 


S 12 545 S12S45 


(3.26) 


I 11 this simple example we have solved a non-trivial integrand just using the building 
blocks and the standard KLT approach. Note that the matrix (3.25) and the vectors 
(3.26) are not simple. Very nicely, one can also use the original KLT approach in 
order to obtain a simpler matrix and vectors [27-29]. For example, let us consider the 
following decomposition 


Isv= v (Wi) i (m KIT ): i 'j(V' i »)|, ( 3 - 27 ) 

oeL,/3ei? 
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where 

L = {( 1, 6(4), 3, 6(5), 2)} = {(14352), (15342)}, (3.28) 

i2 = {(1,4(4),4(5),2,3)} = {(14523), (15423)}, (3.29) 

and 


(UZl)& = m(12345|l, 6(4), 3, 6(5), 2), 

( VIr)o = f dfjL 5 -X- . 

p J (12) (345) (1, /3, 23) 

It is trivial to show that the (^ KLT ) <3 , / § matrix is given by 



and the vectors 


(U?l)& = 

(VI R ) $ = 


(— 

\Sl2 S34 

(— 

\Sl2 S45 


- 1 - 

S12S34 Sl5 s 34 

-1 


Sl 2 S 45 


(3.30) 

(3.31) 


(3.32) 


(3.33) 


so the computation of (3.20) becomes trivial. One can easily check that 

/*» (12345) (12X345) = ? (< 3 - 34 > 

a£L,j3eR 

■Sl5 s 24 
S34 

which agrees with (3.24). 

Finally, consider the following six-point integral 

/^ 6 ( 123456) (12) (34) (56) ‘ 

It turns out that choosing a KLT basis of (n — 3)! = 6 permutations where three labels 
are fixed to some particular locations is not enough in this case. One can show that it 



( 


S 14 S 25 


S 12 S 45 V S 12 S 34 
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is not possible to find such a set so that when multiplied with (12) (34) (56) always gives 
a product of two six-point Parke-Taylor factors. In section 7 we present an explicit set 
of six permutations that do give rise to two Parke-Taylor factors in all six cases and 
hence a building block. 


4 Petersen’s Theorem 


In this section we present the second key result needed for the algorithm in section 
6. Consider a special class of integrands F(a) defined by rational functions with no 
zeroes. This means that F(cr) has 2 n factors cr a b in the denominator and a trivial 
numerator which can be taken to be unity. We denote this kind of integrands as F d (a). 
The required PSL( 2, C) transformation of F d (a) implies that each label a must appear 
exactly in fonr factors. In this section all integrands F satisfy these special properties, 
i.e F{o) = F d (a). 

An integrand with these properties is uniquely determined by a graph (Gp) with n 
vertices constructed by including an (non-oriented) edge connecting a and b for every 
factor of <j a b in the denominator of F d (a). Note that the same factor can appear more 
than once and hence the graph is not simple in general 5 . To illustrate this morphism 
we give a simple example. Let us consider the rational function 


(123)(345)(561)(246)’ 


(4.1) 


Using the previous rules, the graph associated with the function F d (a), which we have 
called Gf, is given by the pair Gf = ( Vf,Ef ), where Vf and Ep are the vertex and 
edge sets, respectively 


V F = { 1,2 ,3,4,5,6}, (4.2) 

E F = {[1; 2], [2; 3], [1; 3], [3; 4], [4; 5], [3; 5], [5; 6], [1; 6], [1; 5], [2; 4], [4; 6], [2; 6]}, 

where [a; b] = [6; a] (non-oriented graph). The line drawing of the graph Gf is 

5 Graphs with multiple edges connecting two vertices are also called multigraphs. 
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1 



(Fig.4.1) Graph G F 

The graph associated with any integrand Fd(c 7 ) has the property that every vertex 
has degree exactly four. These are called 4-regular. In general a graph where each 
vertex has degree k is called fc-regular. 

A result of Petersen dating back to 1891 states that every 4-regular graph G with n 
vertices contains a 2-factor. A 2-factor is a 2-regular subgraph of G with n vertices. Of 
course, after removing all the edges in G from such a 2-factor one is left with another 
2-regular graph with n vertices [30, 31]. 

It should now be clear why this theorem is useful in our construction. Given any 
function Fd(cr) and its associated graph G F , find the two 2-regular subgraphs implied 
by Pertersen’s theorem 1 ’ and denote them by Gp and Gp. Then it is possible to write 

F d = Fi x F «, (4.3) 


where 


F L — 

r d ~ 


n 


£Gp 


O, 


Ve t'U'e 


rpR _ 


fi eS G« 


cr. 


Ve j^e 


(4.4) 


with v e and u e the end vertices of the e edge. Given that Gp and Gp are 2-regular, Ff 
and Fj 1 both transform as half integrands and the generalized KLT construction can 
be used to decompose 


dUnFd 


(4.5) 


as a sum over product of simpler integrals. The concept of “simpler” will be made 
precise in sections 5 and 6. 

6 This decomposition is not unique. 
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For example, in the graph represented by (Fig.4.1) it is simple to see the following 
two 2-factors 



(Fig.4.2) 2-Factor decomposition of Gp. 

and therefore 

rpL _ 7 qiR ___ 

d ~ (123456)’ d ~ (135) (246)' 

Although this example is very simple, finding the decomposition into two 2-factors of 
a general 4-regular graph can be a daunting exercise. In appendix A we provide an 
explicit algorithm for finding two 2-factors Gp and Gp and a sketch of the proof of 
Pertersen’s theorem. 


5 Hamiltonian Decomposition 

In this section we present the third and final ingredient needed for the general algorithm. 
Recall the example presented in section 3.3 

J d ^ 5 (12345)(12)(345)’ 

The KLT procedure worked because we were able to find Parke-Taylor factors such 
that 


(12)(345) x (14523) = (12354) x (12543), (5.2) 

(12)(345) x (15423) = (13542) x (12345). 
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It is clear that by combining Petersen’s theorem with the KLT construction one would 
be able to compute any integral with a trivial numerator if a decomposition of the form 
(5.2) was always possible. This is the main subject of this section and we start by 
introducing some terminology standard in the graph theory literature. 

A 4-regular graph G is said to have a Hamiltonian decomposition if G has two 
edge-disjoint Hamiltonian cycles Hi and H 2 [30]. 

Let us remind the reader what a Hamiltonian cycle is. In any connected graph 
with n vertices one can ask if there is a connected closed path that visits all vertices 
exactly once. Such a closed path, i.e., collection of edges, is called a Hamiltonian cycle. 

We are going to assume that we are given a generic 2-regular graph G R with n ver¬ 
tices, for example when n — 5 we can consider the graph associated to the denominator 
(12) (345). If the 2-regular graph G R is connected then we do not have to continue 
since it is already of the form needed. Assume that the 2-regular graph is made out 
of m disconnected 2-regular graphs, each with the corresponding Parke-Taylor factor, 
M(r 2 ) • • • (r m ). 

We say that an n-point Parke-Taylor factor (a(l)a(2) • ■ ■ a(n)) is compatible with 
a general combination (r!)(r 2 ) ■ ■ ■ (r m ) if the union of both graphs, which is obviously a 
4-regular graph, has a Hamiltonian decomposition. In our more physical terminology, 
(a(l)a(2) • • • a(n)) is compatible with (ri)(r 2 ) ■ ■ ■ (r m ) if 

Or)(r 2 )---(r m ) x (a(l)a(2) ■■■a(n)) = (/3(l)/3(2 ) ■ ■ ■ J3(n)) x (y(l)y(2) • • • y(n)) (5.3) 
for some /3,'y £ S n /Z n . 

How many permutations a £ S n /h n are compatible with a given form (ri)(r 2 ) ■ ■ • (r m ) 
seems to be a complicated question in general. However we have made an extensive 
computer search and have found strong evidence that not only the number is always 
larger than (n — 3)! but it becomes much larger than (n — 3)! as n increases. In the 
next subsection we discuss the results obtained in the computer search as well as a 
conjecture for the precise number when every (r a ) in (rx)(r 2 ) • • • (r m ) contains exactly 
two elements. 

We restrict our search to the cases where all (r a )’s except one have exactly two 
elements (i.e., are bubbles) and the cases where all (r a ) except one have exactly three 
element (i.e., are triangles). 
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5.1 Bubbles and One Polygon 

Let us start with 2-regular graphs of the form 

(12) (34)... (2 m - 1, 2m) (2m + 1, 2m + 2,..., 2m + k). (5.4) 

The graph is then given by m bubbles and one polygon with k sides. The total number 
of vertices is n = 2m + k. 

The results from our computer search are presented in table 5.1 . 


\<6 

X'P 

By 


A 




o 

o 

# 

© 

• 

1 

1 

6 

32 

200 

1,518 

13,160 

126,032 

1,325,682 

2 

16 

144 

1,096 

9,440 

92,712 




3 

480 

6,336 

65,184 

723,360 





4 

28,416 

491,904 









5 

2,588,160 










6 

347,996,160 











Table 5.1. Number of compatible Parke-Taylor factors to the graph obtained by 
combining bubbles (vertical-B) with various polygons (horizontal-Poly). 

Recall that the numbers presented in the table are the total numbers of permuta¬ 
tions that are compatible with the corresponding 2-regular graph. It is easy to see that 
the number in the table are always greater than (n — 3)! = (2m + k — 3)! (except for 
n=4 when it is equal) and that the ratio increases as n gets larger. 

Quite nicely, we have been able to find a sequence of numbers that reproduces the 
case k = 2 for all m tested and this is why we conjecture the following. 

Consider the 2-regular graph with 2s vertices which is made out of s bubbles (here 
s — m+ 1 since the polygon added to the m bubbles in the table is also a bubble). The 
number of compatible Parke-Taylor factors is given by 

2 s-2 (s — l)!A(s) , (5.5) 
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where A( s ) is the number of types of sequential s-swaps moves for the travelling sales¬ 
man problem which has a closed formula presented as sequence A001171 in the OEIS 
webpage [41, 42], 

It is interesting to note that the case of only bubbles is the one with the least 
number of compatible permutations for a given number n of vertices. Using the closed 
formula, it is possible to find the ratio of (5.5) to (n — 3)! = (2s — 3)! as s goes to 
infinity, 


r(s) 


2 s-2 (s — l)!A(s) 7T 

(2s-3)! 2 


(5.6) 


5.2 Triangles and One Polygon 


Finally, we consider the case in which we have several triangles and only one polygon 
with k vertices. Clearly the case of one triangle and k — 2 is already in the table 5.1. 
The results for this new computer search are presented in table 5.2. 


T\ 


A 




o 

o 

# 

• 

1 

6 

42 

264 

1,920 

16,284 

154,098 

1,607,296 

2 

1,404 

14,760 

141,264 





3 

1,111,968 










Table 5.2. Number of compatible Parke-Taylor factors to the graph obtained by 
combining triangles (vertical-T) with various polygons (horizontal-Poly). 

Once again, it is clear that the number of compatible Parke-Taylors is larger than 
(n — 3)! and the ratio grows with the number of vertices. 

6 Main Algorithm 

In this section we present the main algorithm for the decomposition of an integral of a 
general rational function F in terms of the building blocks m(a\/3). 
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6.1 Reduction of the Numerator 


A general function F can have a numerator which can be taken, without loosing of 
generality, to be a monomial in the variables cr ttb . Let us define cross ratios to be the 
basic PSL( 2, C) invariant functions 


_ (Jab&cd 

'f’abcd, = 

®ad&bc 


It is clear that by choosing any pair of Parke-Taylor factors, say 


( 6 . 1 ) 


{( 7 (1) 7 (2)... 7 (n)), (123-■■«)}, (6.2) 

to multiply F, one gets an SL(2, C) invariant function which can be expressed in terms 
of cross ratios, i.e., 


(t(1)t(2) • • -7W)(123- • -n)F = \\r aiM , ClAl . (6.3) 

i =i 

for some value of m > 0. 

Using this we choose as our starting point a representation for F of the form 

F ~ ( 7 (1) 7 (2) ■ ■ ■ 7 ( rt ))(123• • • n) T\ rai ’ bl ' CI ' dl ' (6 ' 4) 

The following is a procedure which might not be the most efficient in particular cases 
but it is general. Let us isolate the m th cross ratio and write 


F = Ft, x F, 


R 


(7(l)7(2)---7 (n)) 


m—1 


X 


(12 ■■■n) 


r a I ,b I ,c I ,di X 


1=1 


^Clm fim ^Cm idrr 
^CLm ,dm ^bm,Cm 


(6.5) 


Lise the KLT procedure to separate the left and right factors (defined by the parenthesis) 
by using a basis of (n — 3)! Parke-Taylor factors of the form l/(c m , d m , a m , a) where a 
is some permutation of the n — 3 left over labels. Applying this procedure one finds 
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that 


/ 


dfAnF' 


m (7(l)7(2)---7('«) 

a,/3 


^(^|/^) / dfd r , 


R 


(c m , d m , 0>rm P'} 


Now we are left with the computation of new integrals in which the factor <J Crnt d m 
in then numerator has been canceled by the same factor arising in the expansion of 
(c m , d m , a m , a). More explicitly we now have to compute 


771—1 


*R 


r ai,bi,ci,d I ) x 


<T, 


firr 


(c m , d m , a m , (3) \ (12 • • • n) j \ <J a m ,d m ( ^b m ,c m (7 dm,a, rn ( ^a rn ,0 i ' ' ' Cr /3n-3,c„ 

( 6 . 6 ) 

Now we can repeat the KLT procedure for each such new integrals. Once again we 
take the left and right factors as those collected in parenthesis in (6.6) but this time 
we use a basis of the form 1 /(a m , b m , d m , a). This gives rise to integrals on the right of 
the form 


dfi r 


® 0 . m-bn 


( cr a m ,d m cr 6 m ,c m cr d m ,a m cr a m ,/3i ' ’ ' ,Cm) ( a mi b m , d m , Cx) 


(6.7) 


which after expanding (a m , b m , d m , a) and canceling the factor cr amtbm in the numerator 
becomes an integral of the special kind discussed in section 4. We leave the computation 
of such integrals to the next subsection. 

On the left side of the KLT formula one has 


^ m— 1 

(a m ,K,d m ,P) (12-.-n) k[ 


( 6 . 8 ) 


But this is identical in structure to the starting point (6.4) but with m — 1 cross ratios. 
Iterating the procedure one ends up reducing the computation of the original integral 



( (7(1)t(2) • • • l( n )) 


X 


771—1 


(12' 


■n 


r a/,6/,C/,d/ X 


1=1 


CT n h CT p r] 

u ‘m i u m '-m ■> Ui n 

CT n rl CT h p 

'- L mi LL m u m Wrr 


(6.9) 


to that of integrals of the special form 



n 


eS G a v e ,u e 


( 6 . 10 ) 
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where G is some 4-regular graph with n vertices. 


6.2 Finding a Compatible KLT Basis 

This is the final step of the algorithm. Consider any integral of the special type as the 
ones found above 



( 6 . 11 ) 


If the integrand can be written as the product of two Parke-Taylor factors then this 
integral is of the form m(a\/3) and it is the end of the procedure. This is the ideal point 
to introduce more of the mathematical terminology corresponding to our physical issue. 
A Parke-Taylor factor of the integrand corresponds to a Hamilton cycle of G. Therefore, 
if G has two edge-disjoint Hamilton cycles we stop. 

Assuming that G does not have two edge-disjoint Hamilton cycles, using Petersen’s 
theorem as described in section 4 we write (6.11) as 



( 6 . 12 ) 


where both Gl and Gr are 2-regular graphs of n- vertices. This means that the left and 


the right factors can be used as the starting points in a KLT decomposition. 

The next step is based on the existence of compatible Parke-Taylor factors as 


discussed in section 5. Recall that a Parke-Taylor factor (a(l)a(2) • • • a(n)) is said to 


be compatible with a 2-regular graph, say Gl, if the union of the two graphs, which 
is a 4-regular graph, admits a Hamiltonian decomposition, i.e., it is the union of two 


edge-disjoint Hamilton cycles (see section 5 for definitions and more details). 


In section 5 we gave computer-based evidence for the fact that for any 2-regular 
graph, the number of compatible Parke-Taylor factors is larger than (n — 3)!. We now 


also assume that from the set of compatible graphs of Gl and of Gr it is possible to 


choose (n — 3)! of each, denoted as C and 1Z such that they are independent in the 
sense defined in section 3.2. For the reader’s convenience we recall that C and 1Z are 
said to be independent of the (n — 3)! x (n — 3)! matrix 


m c ^ n _ { m (a\f3) : a G C, (3 G 7Z} 


(6.13) 


is not singular. 
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Under these assumptions we take C and 77 in order to build a generalized KLT 
relation that expresses 


dp r 


x 


rieeGi a Ve,Ue n e£G R °7e, 


(6.14) 


as a sum over product of 


and 


d/ir 


dfi r 


(IleeGz, cr v e , Ue )(a{l)a(2) ■ ■ ■ a(n)) 


(6.15) 


(6.16) 


' (FUg, ^e,«e)(/^( 1 )/ 3 (2) • • • Pin)) 

where a and f3 define graphs that are Gl and Gr compatible respectively. 

By definition of and Gr compatibility both integrals (6.15) and (6.16) can be 
written as the product of two Parke-Taylor factors, e.g, 


( 11 a ^,u e ){a(l)a(2) ■ ■ ■ a(n)) = (p(l)p(2) • • • p(n))(p'(l)//(2) ■ ■ ■ p'[n)) 


(6.17) 


for some p, p' G S n and therefore 

J dp 

as desired. 


(riesGi a v e ,u e )(a(l)a(2) ■ ■ ■ a{n )) 


= ™(p\p') 


(6.18) 


7 All Six-Point Integrals 

In this section we illustrate the techniques introduced in previous sections in the case 
of six-point integrals. Starting with a general integral 

J dp e F(a) (7.1) 

one can use the reduction procedure explained in the previous section to turn it into 
combinations of the special integrals studied in section 4 which have a constant numer- 
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ator. More explicitly, we are left with the computation of integrals of the form 


d/iQ 


n 


Wab 


l<a<b<6 Gab 


(7.2) 


with w ab > 0, w aa = 0 and J2 b =i "'.a = 4 for all a e {1,2,..., n}. 

As explained in section 4, the matrix u> ab , with = u> ab , is the adjacency matrix 
of a 4-regular graph Gf [30, 31]. Employing Petersen’s theorem any such integrand 
can then be decomposed as 



(7.3) 


where Gp and Gp are two 2-regular graphs that provide one of the possible decompo¬ 
sitions of Gf. 

Each one of the 2-regular graphs Gp and Gp is one of the following options: a 
hexagon; a square and a bubble; two triangles; or three bubbles. 

The last step is to find two sets of six Parke-Taylor factors C and 7Z such that all 
Parke-Taylor factors in C (71) are compatible with Gp (Gp). Of course, the two sets 
C and 7 Z must also be independent. Once this is found, the KLT procedure of section 
3 completes the computation in terms of the basic building blocks m(a\/3). 

At this point is clear that our only remaining task in this section is to provide 
explicit sets of Parke-Taylor factors which are compatible with each one of the possible 
2-regular graph. 

Clearly, there is nothing needed when the 2-regular graph is a hexagon as it is by 
itself a Parke-Taylor factor and therefore compatible with all 5! Parke-Taylor factors 
at six-points. 


7.1 A Bubble and a Square 

Consider the configuration of labels for a bubble and a square given in figure 7.1 
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(Fig.7.1) A bubble and a square geometry. 

Clearly, any other assignments of labels can be obtained from this one by a simple 
relabelling. An independent basis of compatible six Parke-Taylor factors chosen from 
the 32 possible ones is given in figure 7.2 







(Fig.7.2) Compatible cycles with a bubble and a square. 

7.2 Two Triangles 

Consider the configuration of labels for two triangles given in figure 7.3 
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(Fig.7.3) Two triangles geometry. 


An independent basis of compatible six Parke-Taylor factors chosen from the 42 possible 
ones is given in figure 7.4 



(Fig.7.4) Compatible cycles with two triangles. 


7.3 Three Bubbles 

Finally, consider the configuration of labels for three bubbles given in figure 7.5 


i 



2 



3 


(Fig.7.5) Three Bubbles. 

From the table 5.1, we know that the total number of 2-regular compatible graphs 
with three bubble is 16. These 16 Parke-Taylor factors are given in figure 7.6 


29 











(Fig.7.6) Compatible cycles with three Bubbles. 
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In the next section we choose an independent basis of compatible six Parke-Taylor 
factors from the 16 possibilities given in figure 7.6, in order to solve the non-trivial 
example (3.35). 

7.4 Explicit Example 

As an explicit illustration of the full procedure we end this section with the computation 
of the integral we left open at the end of section 3.3 

f 1 1 

4 P - J dpt ( 123456 ) ( 12 )(34)(56) ’ ( 7 ' 4 ) 

with its associated graph 



(Fig.7.7) Graph representation of the integrand T(a) = 2T(ct)Zk(<t), where 
-El(c) = ( 123456 ) ’ = ( 12 )( 34 )( 56 ) 

On the r.h.s the canonical Parke-Taylor (I l{&)) does not have any incompatibility 
with the standard KLT approach, but, for the three bubbles (I r(ct)) the standard KLT 
construction is not enough. The idea is to hnd a set of six cycles (six Parke-Taylor) such 
that the union of each one of them with the three bubbles graph can be decomposed 
in two-disjoint Hamiltonian cycles (i.e. in two Parke-Taylor). 

From the figure 7.6 we choose the following set 

77 = {(145326), (145236), (164523), (153246), (154623), (154236)}. (7.5) 

In order to check that these six elements are compatibles with the three bubbles we, in 
figure 7.8, give an Hamiltonian decomposition 
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1 



(Fig.7.8) Hamiltonian decomposition of (12) (34)(56) (/3), (3 E TZ. The red lines give one 
Hamiltonian cycle and the black lines form the other one (disjoint cycles). 

Therefore, the integrand 

1 

(12) (34) (56) (/3)’ 

with (3 E 1Z, can be written as the product of two six-point Parke-Taylor factors and 
hence their integrals are part of the building blocks. 

So far, we have just found a “right base”, i.e. a linearly independent' TZ set given 
in (7.5). Now, it is necessary to find a “left base”, i.e. a C set which must satisfy the 
following two conditions: 

(1) the matrix defined by the inner product m c,7l (a\/3), with a E C and /3 E TZ, is not 
singular; 

7 Remember, the set is linearly independent using the m(a\(3) inner product. 
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(2) the a's elements, a G C, must be compatibles with the left integrand X l(ct), i.e. 
they admit a Hamiltonian decomposition (see section 5), such as it happened on the 
“right side”, (see Fig.7.8). 

In our example (7.4), the left integrand, X L (cr) = is the canonical Parke-Taylor, 

therefore the second condition is automatically satisfied (product of two Parke-Taylor). 
If we defines the left set (£) with the same elements of the right set, i.e. C = 77, then 
one can show that the matrix m c,7l (a\j3) is not singular. However, note that this is a 
huge matrix since that the diagonal elements are Parke-Taylor squared 8 . We define 
the left set such that its intersection with the right base is disjoint (£ D 77 = 0), in 
order to avoid terms with Parke-Taylor squared in the matrix. For example, choosing 
the linearly independent set 


C = {(126534), (125634), (125643), (124365), (124653), (136524)}, (7.6) 


n 'jp 

one can check that the matrix m a ^ ,qg£ and /3 G 77, given by 


/ £[ 14 : 345 ] 

S35 S26 

0 

0 
0 
0 

V 0 


0 

— (si4 S 25 S 36 ) 1 

0 

0 

0 

0 


0 

0 

£[ 25 : 123 ] 

S13 s 46 

0 

(S 13 S 46 s 123) 1 

0 


0 

0 

0 

(■524 Sl5 ^234) _1 

0 

£[36:234] 

■515 s 24 


0 

0 

(S13 S46 S123) 

£[ 46 : 24 ] 

■535 S135 
£[ 123 : 246 ] 
S13 S46 

0 


0 \ 

£[ 14 : 136 ] 

S25 «36 

1 (S13 S 25 S245) 1 

0 

0 

(^24 ^36 S245) l J 


8 c = n is a good left base if the left integrand is also given by three bubbles, i.e. for a total 
integrand given by the graph 


1 

f 
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1 


1 


(7.7) 


with s ili2 ... j. = (k h + k i2 H-h k tj ) 2 and 

B[ii ■■■i k ■ j i ■ ■ ■ Jm] = 


-1- 

s h--ik s jl--jn 


is not singular. Moreover, the vectors defined as 


{ux c ) a = 

(123456) (a)' “ 6 £ ’ 

(7.8) 

(VX K h = _ 

/ dA6 (12)(34)(56)(/3)’ /?€K ' 

(7.9) 


are given by 

(UI c ) a = 
(Vi*)? = 


B[56:345l / \_l B[34:123l B[12:234l / W 

S12S3 4 ’ -( S 12S34S 56 ) , S12S56 , S34S56 , (Sl2S 56 Sl23) , OJ , 

B[56:345] _/ \_i B[34:456] £[12:234] B[34: 123] B[12:234] \ 

512 534 ’ l S 12 S 34 S 56j ) 512 ^56 ’ 534 S 56 ’ 512 S 56 ’ 534556 J 


Finally, we can write the answer of the (7.4) integral as a rational function of the 
m(a\j3 ) building blocks 


7 sp / d/u, g 


(123456) (12) (34) (56) 


E m% (vz R )„. p.io) 

/3g7\!. 


We have been able to solve a non-trivial integral using the algorithm described 
in this paper. In addition, one can note, such as it was done in the Hamiltonian 
decomposition section, that it is possible to split the geometry of the 4-regular graph 
(initial integrand) and to study them separately as two 2-regular graphs. 


8 Discussions 

In this paper we provided an algorithm for the computation of contour integrals of the 
form 

[ dii n F(a ) ( 8 . 1 ) 
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where F(cr) is a general rational function with the SL( 2, C) transformations of two 
Parke-Taylor factors and therefore can be written, without loss of generality, as 


F(a) 


_-_TT 


( 8 . 2 ) 


where are cross ratios of the positions of four punctures and m is an arbitrary 
positive integer. These integrals appear in many physical applications and hence having 
an algorithm for their computation is important for the study of a variety of theories. 
The contour and measure are defined using the critical points of a Morse function on 
Mo,n [13, 18, 24] 

(p(a,a) = - ^ s ab ln\cr a -a b \ 2 . (8.3) 

l<a<b<n 

This function is a universal part of the “action” that controls correlation functions that 
compute scattering amplitudes in string theory. One of the most pressing issues is to 
find a direct connection with string theory computations. It is well known that in the 
Gross-Mende limit of string amplitudes (which is a certain tensionless limit and thus 
opposite from held theory) correlations functions also localize to the critical points of 
(8.3) [18]. Constructions based on ambitwistor space have been very successful but still 
do not provide a direct link with a limit of string theory [3, 6, 10, 13]. The approach 
of Berkovits that uses an infinite tension limit of string theory in the pure spinor for¬ 
malism is clearly connected to string theory but is not directly connected to the CHY 
formulas for gravity and Yang-Mills [5, 7, 43]. It is reasonable to hope that a direct 
connection of CHY formulas to a limit of string theory may teach us new lessons on 
how string theory is connected to held theory and how the tools developed in this work 
can extend to applications in string theory. 

A simple byproduct of this work is its application to produce a variety of held theo¬ 
retic relations. One of the simplest examples is a relation that expresses double-trace 
amplitudes of gluons in the Einstein-Yang-Mills theory discussed in [11, 12] to single 
trace amplitudes and scalar amplitudes [1, 2], To see the relation consider the CHY for¬ 
mula for a double-trace amplitude with gluons 1, 2,..., m in the hrst trace and gluons 
m + 1, m + 2,..., n in the second, 


A^ 2 \l, 2 ... , m : m + 1, m + 2,..., n) 


dfi r 


Sl2---r 


(12 • • • m ) (m + 1 m + 2 • • • n ) 


Pf'T. 
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The graph of the integrand is simply the product of two polygons and hence one can 
find a basis of (n — 3)! permutations (Parke-Taylor factors) that are compatible with 
that graph. Recall that by being compatible we mean that the union of the graph 
associated to a permutation and the two polygons admits an edge-disjoint Hamiltonian 
decomposition. In other words, it can be expressed as the union of two Parke-Taylor 
factors. Denoting such as a basis as a left basis C and choosing any other convenient 
basis as a right basis 1Z one has 


A (2) (l,2... ,m : m+l,m+2,... ,n) = s 12 ... m m(a'\a") (m cln ) al0 A (1 \(3), (8.4) 

a£C,0£TZ 


with 


m(oi'\a") = / 


(a) (12 • • • m) (m + 1 m + 2 • • • n ) 


/ dfi n . ' 


a') (a' 1 


(8.5) 


and where a! and a" are the two edge-disjoint Hamilton cycles that decompose 


(a) (12 • • • m) (m + 1 m + 2 • • • n ). 

It would be interesting to explicitly construct the bases C and 1Z such that the relation 
(8.4) takes its simplest possible form. Also interesting is to find a possible string 
theoretic origin for this relation. It is clear that there are plenty of relations such as 
(8.4) which connect very different kinds of objects. In fact, these new kind of relations 
greatly extend the large class already found using KLT in [12] which linked theories 
such as the special Galileon with the U(N) nonlinear cr model. The algorithm presented 
in section 6 is completely general but it is not the most efficient one in particular cases. 
It would be interesting to select certain families of integrals that appear in particular 
theories and refine the algorithm to make it the most efficient possible. One way to 
improve the efficiency is by selecting basis of permutations C and 1Z such that the 
matrix is as sparse as possible. Finally, finding direct mathematical application 

of contour integrals over A4o,n could lead to yet another link between the elegant math 
of Riemann surfaces and that of graph theory. One closely related link was established 
by using Strebel differentials [44] which provide a connection between the decorated 
moduli space, -Mo,n x (M + ) n and ribbon graphs [45]. It would be fascinating to find 
connections with this or other constructions. 
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A Algorithm and Petersen’s Theorem 

Before stating the Petersen’s Theorem it is useful to give an algorithm to find a de¬ 
composition of any 4-regular graph in two edge-disjoint 2-factors, in order to elucidate 
the Petersen’s Theorem. 

A.l Algorithm 

The main aim of this algorithm is to show in a simple way how to obtain two edge- 
disjoint 2-regular spanning 9 graphs from one 4-regular graph, i.e. 

G f = G l f u 

where the “U” symbol means the union of the edge sets 10 . 

The algorithm is divided in four fundamental steps. In order to better understand 
it we give an example in each step using the graph drew in the (Fig.4.1). 

• Eulerian Path (EP) 

The first step of the algorithm is to find an Eulerian path. An Eulerian path 
is a trail in a graph which visits every edge exactly once. Since we are working 
with 2/c-factor undirected graphs, k G N+, then the existence of an Eulerian path 

9 A Spanning graph is a subgraph such that its vertex set is the same as the original graph. 

10 The vertex sets are the same for the two 2-regular graphs and the 4-regular graph, i.e. Vp = Vp = 
Vp. 
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is guaranteed and, in addition, every Eulcrian path is closed, i.e. it is a cycle 


[30, 31]. 


An Eulerian orientation of an undirected graph Gf is an assignment of a direction 
to each edge of Gp such that, at each vertex a*, the indegree of a* (number of 
incoming edges) equals the outdegree of a, (number of outgoing edges). 

There are several algorithms to get an Eulerian path but we do not discuss them 
in this paper [30, 31]. In order to unify the notation we denote the cycle defined 
by the trail {[ap a 2 ], [ap a 2 ],..., [a„_i; a n \, [a n ; oq]} as 


(fllCZ 2 .. -®n) — { [fll, Ct 2 ] j [®2; O3] j ■■ - j [®n.— 1j j [®m ®l] 


where the orientation of the path is given by the ordering of the edges in the set. 

With a view to understand better the Eulerian path concept we have found two 
of them on the example given in the (Fig.4.1) 



(Fig.A.l) Eulerian paths. 


Note that on each vertex there are two incoming and two outgoing edges, such 
as it was said above. 

• Bipartite Graph 

After finding an Eulerian path, the second step is to construct a bipartite graph 
from this EP. Thus, we must define what is a bipartite graph. 

Bipartite Graph 

A bipartite graph G is a graph whose vertex-set V can be partitioned into two 
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subsets U and W, such that each edge of G has one endpoint in U and one 
endpoint in W. The pair U, W is called a (vertex) bipartition of G, and U and 
W are called the bipartition subsets [30, 31]. 

Using the EP found in the Erst step one can associate a bipartite graph G' F to 
the original graph Gf- The idea is simple, since on each vertex there are two 
incoming and two outgoing edges then one can split each vertex in two, one of 
them contains the incoming edges (we call them the blue vertices) and the other 
one contains the outgoing edges (we call them the yellow vertices). This procedure 
generates a bipartite graph G' F , where the bipartition is given by the blue and 
yellow vertices, i.e. 


U' F = {Blue Vertices}, W F = {Yellow Vertices}, 

and therefore each edge of G' F has one endpoint in U' F and one endpoint in W F . 

Using this procedure on the Eulerian paths given in the Figure (Fig.A.l) one 
obtains the following 

Bipartitie Graph 




(Fig.A.2) Bipartite graphs 

where we have kept the orientation of the edges on the bipartite graph, but this 
is not necessary. Note that the two bipartite graphs in (Fig.A.2) are 2-regular. 
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• 1-Factor (Perfect Matching) 

The third step is to identify two edge-disjoint (disjoint sets of edges) 1-Factors 


in the bipartite graph G' F . By definition, a 1-Factor in a graph G is a 1-regular 
spanning subgraph, i.e. there is a set of edges without common vertices and 
connects all vertices of the graph G (perfect matching). Moreover, since we have 
found a bipartite 2-regular graph in the step 2, then there is a theorem which 
implies the existence of a perfect matching in it [30, 31]: 

Theorem A.l 

Every r-regular bipartite graph G with r > 0 is 1-factorable, i.e. there are r 
edges-independent 1-Factors 11 . 

This powerful theorem means that the bipartite 2-regular graphs G' F obtained 
in the step 2 have two edge-disjoint 1-Factors. These 1-Factors can be obtained 
easily. First of all, one should detect all cycles in the bipartite 2-regular graph G' F 
and after that, one should remove alternate edges in each cycles. This procedure 
generates one 1-factor and the other one is given by the edges removed. 

In the bipartite graph given by the (Fig.A.2)(a) we have found two edge-disjoint 


1-Factors 


1-Factor 



MS ©— 

®—K?) „ 

© ® 


(a) 


© © 



(Fig.A.3a) Two edge-disjoint 1-Factor graphs 


n This theorem is not proved in this paper. 
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and in the bipartite graph (Fig.A.2)(b) we choose the following two edge-disjoint 
1-Factors 

1-Factor 







(Fig.A.3b) Two edge-disjoint 1-Factor graphs. 


• Identifying Vertices 

The last step is very simple. The idea is just to identify the blue and yellow 
vertex with the same label in each 1-Factor graph, i.e. 



(Fig. A.4) Identifying Vertices. 

and so one obtains the decomposition of the 4-regular graph Gp in two edge- 
disjoint 2-regular graphs ( Gp,Gp ). 

In the figures (Fig.A.5a) and (Fig.A.5b) we show graphically the whole algorithm 
using the two Eulerian paths found in the example of the step 1. 
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Eulerian Path 



Bipartite graph 


1-Factor 





(Fig.A.5a) 
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Parke-Taylor (153246) 


(126543) 


(Fig.A.5b) 


Note that we have obtained two different decompositions. The first one (Fig.A.5a) 
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is one Parke-Taylor factor and two triangles 


Fd (123) (345) (561) (246) Fd Fd ' Fd (135) (246)’ Fd (123456)’ 
and the second one (FigA.5b) is two Parke-Taylor factors 

rp _ ___ _ rpL piR jpL _ _ \ _ rpR _ ___ 

d ~ (123) (345) (561) (246) “ d d ’ d ~ (153246)’ d ~ (126543)' 

A.2 Petersen’s Theorem 

Now, we are able to state the Petersen’s Theorem [30, 31], which was sketching previ¬ 
ously with the algorithm. 

Petersen’s Theorem 

Every 2/e-regular graph G, with k 6 N + , is 2-factorable. 

In the particular case when k — 2 one has that every 4-regular graph is two fac¬ 
torable, i.e., there are two edge-disjoint 2-regular graphs. 
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